Collective fluorescence and decoherence of a few nearly identical quantum dots 
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We study the collective interaction of excitons in closely spaced artificial molecules and arrays 
of nearly identical quantum dots with the electromagnetic modes. We discuss how collective flu- 
orescence builds up in the presence of a small mismatch of the transition energy. We show that 
a superradiant state of a single exciton in a molecule of two dots with realistic energy mismatch 
undergoes a two-rate decay. We analyze also the stability of subdecoherent states for non-identical 
systems. 



I. INTRODUCTION 

Confinement of carriers in semiconductor quantum 
dots (QDs) leads to spectrally isolated states which may 
be optically controlled at a high level of coherence^. 
A single QD offers at most two degrees of freedom (a 
biexciton) which may be coherently manipulated by opti- 
cal fields in various ways 3 ""^, allowing one to demonstrate 
the simplest non-trivial quantum logical operations. In 
order to overcome this two-qubit limitation one needs 
to develop manufacturing methods and control schemes 
for arrays of two and more QDs. Conditional control of 
such systems, indispensable in both classical and quan- 
tum computing schemes, requires interaction between the 
QDs in the array^, which may be provided, e.g., by 
dipole interaction between confined excitons^. There- 
fore, much experimental effort has been devoted to the 
investigation of the coupling between QDs and its sig- 
natures in the optical response and correlation statistics 
of quantum dot molecules (QDMs) built of two coupled 
QDs^ ' 10 ' 11 ' 12 . It turns out that dephasing of excitons in 
QDMs differs considerably from that of individual QDs 13 . 

Even without interaction, optical properties of QDMs 
and QD arrays may be strongly modified due to collec- 
tive coupling of sufficiently closely spaced QDs to the 
electromagnetic (EM) field. These collective effects have 
been extensively studied for atomic system a 14 ' 15 ! 16 ] 17 ? 18 
where they manifest themselves by superradiant emis- 
sion, i.e., an outburst of radiation from the excited sam- 
ple, markedly different from any exponential decay^. 
On the other hand, the collective interaction leads 
to the appearance of subradiant states for which the 
probability amplitudes for photon emission interfere de- 
structively, leading to decoupling from the EM reser- 
voir and to infinite lifetime. It has been proposed 
to use these states for noiseless encoding of quantum 
information 2 ^. For the analogous problem of coupling 
to phonon modes of a semiconductor, "subdecoherent 
states" of QD arrays have been suggested as a possible 
noiseless implementation 21 . Compared to atomic sam- 
ples, QDs may be easier to arrange in a regular array 
but the perfectly identical transition energy characteris- 
tic of natural atoms is extremely hard to reach for these 
artificial systems. 

The purpose of this paper is twofold. First, we deal 
with the general, theoretical problem of the stability 



of collective interaction, including noiseless encoding, 
against variations of the transition energies. Second, we 
look for clear manifestations of interaction between the 
QDs in the experimentally observable optical properties 
of the QD arrays, depending on the difference (mismatch) 
of the transition energies of individual QDs. 

Thus, we study the interaction between small, slightly 
inhomogeneous arrays of QDs and their EM environment. 
The system evolution is described within the Weisskopf- 
Wigner approach 2 ^. We show how the coherent interac- 
tion is destroyed by growing inhomogeneity of the tran- 
sition energies in the regime where the latter is compara- 
ble to decay rate. As important examples, we discuss the 
buildup of superradiant emission from a few QDs and the 
decay of "subdecoherent" states built on non-identical 
QDs. We show that interaction between dots in a regu- 
lar array may, to some extent, stabilize the coherence, in 
contrast to randomly distributed atomic systems, where 
it leads to dephasingi^. We discuss also how the inter- 
play of the transition energy mismatch and interaction 
strength determines the time evolution of the lumines- 
cence of a more realistic QDM. We show that the decay 
of luminescence varies from exponential with a single QD 
rate (for weakly interacting dots), through nonexponen- 
tial (for interaction compared to energy difference), again 
to exponential with doubled rate (when interaction en- 
ergy prevails). 

The paper is organized as follows. In Section [IT] we 
define the model used to describe the system. Section 
HTT1 describes the collective flourescence and stability of 
quantum states in a molecule built of two dots. Next, 
in Section IIVI we extend the discussion to arrays of four 
dots. The final discussion and conclusions are contained 
in Section [V] 



II. THE SYSTEM 

We consider an array of QDs located at points r; . We 
assume that each dot may either be empty or contain 
one ground state exciton of fixed polarization with an 
interband transition energy Ei, hence can be described 
as a two-level system. The dots interact with transverse 
EM modes with frequencies LOf. — ck, where k is the 
wavenumber and c is the speed of light. We will describe 
the system in the interaction picture with respect to the 
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Hamiltonian H = E^j n>j + J2k,\ w k&k,A&k>A ( fi = 1)j 
where &fc,Aj ^j- \ are photon creation and annihilation op- 
erators (A labels polarizations), E is the average of the 
energies Ej , and hj is the occupation operator for the j th 
dot, fij = a\ <j_ , where CT w = ,#'±10$" and 

^;c,y &re 

Pauli matrices acting on the jth two-level system. The 
Hamiltonian of the system is then H = Hx + Hi. The 
first component describes the excitons, 

3 ¥j 

where Aj = Ej — E are the energy deviations from the 
average and Vij are Forster couplings between the QDs 34 , 

y. 1 U ^ d -^ \ rr-r- ri 

where ti is the interband dipole moment (for simplicity 
equal for all dots), eq is the vacuum dielectric constant, 
and e r is the relative dielectric constant of the semicon- 
ductor. For self-assembled dots, typical values for this 
coupling range from /xeV for the distance ry of order 
of 100 nm to meV for closely stacked dots separated by 
~ 10 n m 12 i 23 i 24 . Another contribution to the coupling 
may come from the polariton effect (coupling to trans- 
verse field) 25 . 

The second term in the Hamiltonian accounts for the 
interaction with the EM modes in the dipole approxima- 
tion and rotating wave approximation (RWA) 

^ = E- W E-9^ fc - £) feL + H.c, (2) 

l k,X 

with gk\ = id ■ e\(k) ^ ujk / {2eqE v v) , where e\(k) are unit 
polarization vectors and v is the normalization volume 
for EM modes 3 - 5 -. The QDs are placed at distances much 
smaller than the relevant photon wavelength so that the 
spatial dependence of the EM field may be neglected (the 
Dicke limit). For wide-gap semiconductors with E ~ 1 
eV, zero-temperature approximation may be used for any 
reasonable temperature. 

In the present discussion, we disregard the coupling of 
the carriers with phonons. Let us note that the quan- 
tum confinement of excitons leads to a separation of at 
least a few meV between the ground exciton state in- 
volved in our analysis and the lowest excited state in a 
single dot. Therefore, no real phonon-induced transitions 
may take place in a single dot as long as the temperature 
is low enough. It has been shown that the combina- 
tion of dipole interaction and phonon coupling may lead 
to phonon-assisted Coulomb transfer between the dots 2 ^, 
which might be responsible for the uni-directional trans- 
fer observed in the experiment 1 - 2 *. However, the estimated 
rate reaches its maximum of ~ 2 ns for the energy sep- 
aration of a few meV and decreases considerably away 
from this point^i. Therefore, we neglect this effect in 
the present considerations. For extremely closely spaced 



dots, with strongly overlapping carrier wave functions, 
phonon-assisted tunneling 27 processes might also take 
place on time scales comparable to those characteristic 
of the radiative decay. Similar to the phonon-assisted 
Coulomb transfer, such processes would lead to thermal- 
ization of the state of a QDM or QD array which, in 
general, might suppress the dynamics described in the 
following sections. 

Another phonon effect on the exciton state is pure 
dephasing 2 ^' 29 . In QDMs, like in individual QDs, such 
processes affect only the first few picoseconds of the op- 
tical response of a QDM 13 , while our present discussion 
is focused on the radiative decay that develops at much 
longer times. Due to this separation of time scales, the 
evolution related to the radiative processes may be dis- 
cussed separately from this pure dephasing effect. If the 
system state is prepared by an ultrafast pulse, the initial 
phonon dynamics may result in a certain reduction of the 
coherence of the initial state, as we qualitatively discuss 
in the following sections. 



III. QUANTUM DOT MOLECULES (2 QDS) 

We will start our discussion with quantum dot 
molecules composed of two QDs. In the present sec- 
tion, we will first discuss the decay of sub- and super- 
radiant single exciton states in terms of the formal quan- 
tum fidelity with respect to the unperturbed state and 
in terms of the experimentally measurable exciton occu- 
pation. Then, we proceed to the decay of the biexciton 
state which will be studied again in terms of fidelity and 
in terms of the measurable photon emission rate. 



A. Single-exciton states 

The RWA Hamiltonian conserves the number of exci- 
tations (excitons plus photons). Let us first consider the 
initial subradiant state \tp(0)) = (|01) — \10))/y/2, where 
the two-digit kets denote the occupations of the respec- 
tive dots. Since there is only one excitation in this state 
it may, in general, evolve into 

M*)> = c oi (*)|01) + cio (t) \10) + J2 c ookx (*)|00, kX) , 

k,X 

where the last ket denotes the state with no excitons and 
with one photon in the mode (k, A). The Schrodingcr 
equation leads to the system of equations for the coeffi- 
cients 

ic 01 = Ac m + Vc w +Y.9*kx^okxe l(E - u ' k) \ (3a) 

fc,A 

ic w = -Acw + Vcai+Y^aixCookxe^-^tSb) 

k,X 

icookx = 9kx(c i +c w )e-< E - u ^*, (3c) 
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FIG. 1: (a) The fidelity for a subradiant state for V = 
and A = 0.8 ^eV (solid), 2 ^eV (dashed) and 20 ^eV (dot- 
ted), (b) The fidelity for V / 0, A = 2 peV. Inset in (a) 
shows the fidelity of a superradiant state for different dots 
without interaction (A = 4 peV, solid) and for interacting 
dots (A = 2.68 peV, V = 2.97 peV, dashed), compared to 
an exponential decay with the rates V and 2r (dotted and 
dash-dotted, respectively) . 



where V = Viz and A = A2 = — Ai. Following the 
standard Weisskopf-Wigner procedure^ we formally in- 
tegrate Eq. (f3"c)) and substitute to Eq. (j3"al b). which yields 

coi,io = T^Acoi - iVcio 

dsTZ(s)(c Q i(t - s) + c w (t - s)), 



where TZ(s) — J2k a \9kx\ 2 e t ^ E ~ UJkS>t is the memory func- 
tion of the photon reservoir. As the latter decays ex- 
tremely quickly compared to the timescales of the evolu- 
tion of 001,10 one can perform the usual Markov approxi- 
mation and neglect the s dependence under the integral. 
Using the fact that 



Re / dsK(s) 



E 3 \d\ 2 

6lTEoC 3 



r 

2 : 



where T is the spontaneous decay rate (throughout the 
paper we set 1/r = 1 ns), neglecting the imaginary part 
(i.e., assuming that the Lamb shift and other radiative 
corrections are included in the energies), and defining 
c = (cio,c i) T one gets 



c = Ac, A = iAa z 



iV ■ 



(4) 



where a x ,z are Pauli matrices and I is the unit matrix. 
The reduced density matrix for the charge subsystem 
may now be easily constructed as poi,oi = |coi| 2 , P1040 = 

I cio 1 2 , An, 10 = Pio,oi — c oi c iO: Poo, 00 = 1 — Poi.oi - Pio,io> 
with all the other elements equal to 0. 

In order to test the stability of the ideally subradiant 
state in the case of non-identical dots we denote by \ip(t)) 
the pure state evolving from \ip(0)) in the absence of the 
EM reservoir (r = 0) and define the fidelity of the actual 
state p by F = \J (%\){t)\p\%l){t)) . In Fig. QJa) we show 
the result for a few values of the energy difference A 
in the limit of vanishing Forster coupling between the 
dots (i.e., for sufficiently distant dots). It is clear that 




FIG. 2: The exciton occupation for sub- (a) and superradiant 
(b) states for A = 1 meV. The inset in (b) shows the values 
of the occupation decay rates. 



the state maintains its subradiant (stable) character until 
t ~ 7r/(2A) but then it enters a superradiant phase and 
the fidelity rapidly decays below the value corresponding 
to an exponential (uncorrelated) decay. Depending on 
the value of A, a certain number of oscillations around 
this uncorrelated decay rate may be observed. In the 
limit of large A these oscillations become very fast, their 
amplitude decreases, and the decay closely follows that of 
uncorrelated systems, as expected for systems with large 
energy difference and therefore interacting with disjoint 
frequency ranges of the photon reservoir. It is clear that 
observing collective effects for such non-interacting dots 
requires transition energies identical up to several peV. 

If the QDs are close enough, the Forster interaction be- 
comes effective. Since the sub- and superradiant states 
are eigenstates of the Forster Hamiltonian separated by 
an energy 2V , the transition from the initially subradiant 
state to the superradiant state is suppressed if the magni- 
tude of the Forster coupling exceeds the energy difference 
A. This is shown in Fig. [TJb). It is clear that the de- 
cay rate is reduced when V ~ A and the subradiance is 
recovered for V ^S> A. Note that, apart from the triv- 
ial limiting cases, the decay is markedly non-exponential 
and its modulation yields information on the origin of the 
energy level splitting in the system. Indeed, the decay of 
the superradiant state |-0(O)) = (|01) + \10))/V2 shown in 
the inset to Fig. QJa) i s clearly different for two systems 
with the same energy splitting, depending on whether 
the splitting originates from the difference between the 
dots or from the interaction. 

The signatures of collective interaction with the elec- 
tromagnetic field may also be found in the evolution of 
measurable quantities. As an example, let us consider the 
average number of excitons in the QDM. In the present 
state of the art of QD manufacturing, the differences be- 
tween the transition energies of the two dots are rather in 
the meV than in the /j,eV range discussed in the previous 
case. Therefore, let us consider the evolution of exciton 
occupations for the initial states |-0(O)) = (|01)±|10))/v / 2 
for a QDM with A = 1 meV. The results are shown in 
Fig. [2] As can be seen, in this case the decay of the oc- 
cupation shows no oscillations. For V <C A, both states 
show simple exponential decay with the rate T. In the 
opposite limit, V 3> A, the subradiant state becomes 
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stable while the superradiant state decays exponentially 
with a twice larger rate. In the intermediate range of 
parameters, the decay is not exponential. 

The superradiant state is particularly relevant for opti- 
cal experiments since such a bright combination of singlc- 
exciton states is excited by ultrafast optical pulses from 
the ground state. Thus, the curves plotted in Fig. [5] di- 
rectly correspond to the decay of population after an op- 
tical excitation: 36 Therefore, let us study this decay in 
more detail. 

The system evolution is governed by the matrix A 
[Eq. Q] which, for r <C A, is nearly purely anti- 
Hermitian. Therefore, its eigenvectors are nearly or- 
thogonal to each other. Up to corrections of order of 
r/A <C 1, they can be written in the form 



u+ 



COS If 

— sin ip 



where 



sm tp = 



1 

71 



i - 



A 



sm<£> 
cos if 



1/2 



VA 2 + V 2 
The corresponding eigenvalues are 

r 



a± = -- ± v/-A 2 + C^ + r/2) 2 . 

The solution of Eq. ^ for the superradiant initial con- 
dition is 



c(t) = sin(ip + 7r/4)M_e + cos(</? + ir/4)u + e 
The number of excitons therefore evolves as 



x+t 



n(t) = \c(t)f 



sin 2 (^ + 7r/4) e 2RcA -* 
+ cos 2 (^ + 7r/4)e 2RcA+t . 



Due to the almost perfect orthogonality of the eigenvec- 
tors u± the interference term vanishes and the occupa- 
tion decay is a combination of two exponentials with dif- 
ferent rates, as shown in Fig. (Hb). In the inset to this 
figure we show the values of the two decay constants as 
a function of V for A = 1 meV. 

If the initial sub- or superradiant state of a QDM 
is prepared by an ultrafast optical pulse it will partly 
loose its coherence within a few picoseconds of the sys- 
tem evolution due to phonon-induced pure dephasing. 
The detailed dynamics of this dephasing process differs 
from that of a single QD and depends on the system 
geometry 31 . Nonetheless, its essential effect is to per- 
turb the superposition state towards a mixture of two 
states, each of which undergoes the usual exponential 
decay. Therefore, one may expect a decrease of the am- 
plitude of the oscillations in Figs. [1] and [3] and a shift of 
the decay curves in Fig. [2] towards the monoexponential 
decay with the usual decay rate. For special values of 
the energy mismatch, the results may also be modified 
by the phonon-assisted Coulomb transfer 26 . 



B. Biexciton state 

Next, let us consider the case of the same two QDs, 
but initially excited to the |11) state. This state can be 
experimentally prepared in various ways 3 -^. The general 
form of the state is now 



|V(t)>= Cll (t)|ll> 



fc,A 



coi fe A(«)|01feA) 



fc,A 



ciofcA(*)|10fcA) 



cook\ q \>(t)\00,k\q\'), 



and the amplitudes evolve according to the equations 

"•u = ? ]9k\ ( c oifcA + c wkx )e l{E ~ UJh) \ (5a) 

,-i(E-u) k )t fefo\ 



fc : A 



icoikx = Acqi + Vc 10 + gk\cne~ 



g q \>cook\ q \>e 

q,X> 



iciokx = -Acio + Vc 01 + gk\c\\e 



-%{E-u> h )t 



(5c) 



+ Y,9 q yCookx g ye^ E -^ t , 

icookXqX' = g q X'(coikx+ciokx)e~ l(E ~ UJ '' )t . (5d) 

As previously, we formally integrate Eq. (|5d[) . insert it 
into Eqs. (|5bl5c[) , and use the short memory assumption. 
This yields the equation for c fcA = (ciofcA, c ifcA) T , 



CfeA 



igkx / dse A ^ Cll (s)e-^ E ^ s b, (6) 



where b — (1, 1) T and A is defined in Eq. (|4]). Substitut- 
ing this in turn into Eq. (|5ap we find 



c n = - / ds^(s)b T e j4s bcn(i - s). 
Jo 

Since the elements of A are of order of /xeV or meV, the 
matrix exponent is slowly varying on the timescales of 
reservoir memory, as is Cn(t), and both may be taken 
at s = 0, which leads to the decay equation in the usual 
form cn = — Ten or, for the corresponding element of 
the reduced density matrix, pii,ii = — 2rpn.n. 

The evolution equations for the other elements of the 
density matrix may be found by writing, for instance, 
p i,oi = 2R- e Z)fc,A( c iofcACiofcA), substituting ciofcA from 
Eq. and using once more the short memory approxi- 
mation. Performing this procedure for all the elements of 
p in the single-exciton sector, one arrives at the equations 

Ai = -2r/u, 

Ao = -T(Rcp + f 10 -f 11 ) + 2VImp, 
Ai = -r(Rep + /oi-/ii)-2UImp, 

p = -2iAp-T(p+ /n) 

+iV(foi - fio), 
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FIG. 3: Photon emission rate for a superradiant state for FIG. 5: Photon emission rate for a superradiant state of 4 
V = (a) and for V ^ 0, A = 2 /ieV (b). QDs for V = (a) and for V / (b). 
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FIG. 4: The fidelity for two subradiant states IV'a) (a) and 
(b) of 4 QDs for various combinations of parameters. 
The parameter values given in (b) are valid for both figures. 



where we denoted fij = pijjj, l,j — 0, 1, and p — poi.io- 
The photon emission rate 7 = — (2/n + foi + /10) for 
the initial state |11) is plotted in Fig. [3l In the case of 
V = [Fig.[3ja)] we see that the photon emission loses its 
superradiant behavior for growing energy mismatch be- 
tween the dots, tending to the usual exponential decay for 
large A. Like in the previous case, removing the degen- 
eracy between the sub- and superradiant single-exciton 
states by including the Forster coupling stabilizes the col- 
lective fluorescence [the dotted line in Fig.[3Jb) coincides 
with the A = line in Fig. ISJa)]. 



IV. QUANTUM DOT ARRAYS (4 QDS) 

In this section, we study arrays of four QDs in a very 
special, regular arrangement. The resulting symmetry of 
the Forster term leads to symmetric eigenstates and, as 
we show below, to the stabilization of collective effects. 

In general, the Weisskopf-Wigner equations lead to the 
Lindblad equation for the evolution of the reduced den- 
sity matrix of the charge subsystem 



P : 



■i[H x ,p]+£\p], 



(7) 



with 



S_pE H 



where E± = J2j a ±*- We now use Eq. ([7]) to study the 
evolution of four QDs forming a square array in the xy 



plane. The energy deviations of individual dots are now 

A, = aiA, where J2i a i = ancl Si a i — 1> so tnat A is 
the mean square variation of the transition energies. The 
details of the system evolution depend on the particular 
choice of on but the general behavior is only governed 
by the interplay of A and V (unless some particularly 
symmetric choice is made). We arbitrarily fix a\ = 0, 
Q'2 = —0.8, as = 0.27, an = 0.54 and use the mean square 
variation A as a parameter. The Forster interaction is 
parameterized by its magnitude V, with V12 = V23 = 
V34 = V41 = V and V13 = F 24 = 2- 3 / 2 y (the dots are 
numbered clock- wise). 

First, let us choose the subradiant initial states 
|^ a (0)) = (|1001) - |0101) + |0110) - |1010»/2 and 
|^b(0)) = (|1001) - |0011) + |0110) - |1100»/2, which 
span the subspace of logical qubit states that may be 
used for noiseless encoding of quantum information on 
four physical qubits^. Obviously, for non-identical dots 
the phases in these superpositions will rotate and the 
state will be driven out of the initial noiseless subspace 
which leads to a decrease of fidelity, as shown in Fig. 2] 
Out of the two states, only 1^(0)) is a non-degenerate 
eigenstate of the Frohlich interaction for the square ar- 
ray. As a result, as can be seen in Fig. |H only this state 
is fully stabilized by the Forster interaction for V ^> A 
(the lines for V = 100 fieV in Fig. 0] are very close to 
the asymptotic case of V — > 00). Since the other state 
IV'a(O)) is never completely stable the entire "noiseless 
subspace" of logical states^ remains stable only for an 
extremely homogeneous array of QDs. 

Finally, let us study the photon emission rate from a 
superradiant state of four excited QDs, \ip(0)) — |HH) 
(Fig. [5]). Now, a clear superradiant peak of photon emis- 
sion develops for identical dots but vanishes as the dots 
become different. Again, interaction between the dots 
in a regular array stabilizes the collective emission. It 
is interesting to note that the superradiant emission is 
close to ideal already for V ~ A, while the subradiant 
states are stabilized only when the interaction exceeds 
the energy difference by an order of magnitude. 

It should be stressed that the stabilization effect re- 
sults from the special, highly symmetric arrangement 
of the QDs. It should be contrasted with the dephas- 
ing induced by analogous interactions in the randomly 
distributed atomic samples 15 . Likewise, in an irregular 
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ensemble of QDs obtained by spontaneous self-assembly 
no stabilization effect should be expected. However, re- 
cent progress in the pre-patterned and strain-engineered 
growth of QDs 32,33 shows great promise for the manufac- 
turing of QD arrays with a desired geometry. 

V. CONCLUSIONS 

We have shown that collective interaction of carries in 
QDs with their EM environment is extremely sensitive 
to the homogeneity of the QD array. Already for the en- 
ergy mismatch of order of fj,eV the sub- and superradiant 
behavior of physical quantities is replaced by their oscilla- 
tion around the average exponential decay. Thus, observ- 
ing collective fluorescence effects in an ensemble of non- 
interacting QDs seems highly unlikely. Likewise, imple- 
menting the noiseless encoding schemes requires the level 
of homogeneity much beyond the reach of the present 



technology. 

The destructive effect of inhomogeneity can be, to 
some extent, overcome by excitation-transfer coupling 
(Forster or tunneling) between the dots placed in a reg- 
ular array. This can stabilize the subradiance of a state 
of 2 QDs and the superradiant emission from 4 QDs but 
(for a square alignment) still cannot assure stability of 
the entire noiseless subspace implemented on 4 QDs. 

When the energy mismatch between the dots is of 
order of meV, like in the currently fabricated artificial 
molecules of two QDs, the oscillations disappear and one 
observes a decay of the excitation (exciton occupation 
number) composed of two exponentials. For such a re- 
alistic energy mismatch, the two decay rates for non- 
interacting dots are practically equal to the free decay 
rate T. However, with growing interaction strength they 
approach 2T (superradiant component) and (subradi- 
ant). Thus, their values carry information on the origin 
of the energy splitting (interaction vs. energy mismatch). 
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